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ABSTRACT 

In this paper we address the problem of joint transmitter and 
receiver optimization in a wireless multiuser MIMO-OFDM 
environment. We consider a centralized design strategy as- 
suming that Channel State Information (CSI) is available. 
The optimization goal is the niinimizatioD of the total trans- 
mitted power subject to different constraints, such as the 
Quality of Service (QoS) for each user in terms of the mean 
Bit Error Rate (BER) or the maximum power for each trans- 
mitter. In this design, the physical resources are used in an 
optimum way and the Multiple Access Interference (MAI) is 
. jniinimized. The proposed technique is a variant of the Sim- 
JJilated Annealing (SA) algorithm, which is able to find the 
global optimum. Besides, it can handle all kind of constraints 
by generalizing the definition of the cost function. Finally, 
we present simulation results and compare this technique 
with a Lagrange-gradient based solution. 

1 INTRODUCTION 

In last years, special attention has been given to the joint use 
of smart antennas at both the transmitter and receiver, con- 
figuring a MultH Input-Multi-Output (MIMO) channel. The 
new standards such as [1], tend to specify very high fre- 
quencies for the front-end and Radio-Frequency (RF) chains, 
permitting the use of multiple antennas at both the Access 
Points or Base Stations (BSs), and the Portable or Mobile 
Terminals (MTs). AdditionaUy, the Orthogonal Frequency 
Division Multiplexing (OFDM) modulation has been pro- 
posed for several broadcast ins and communications systems 
such as the European WLAN Hiperlan/2 [1]. 

There is also an increasing demand for higher capacity 
or number of users, and Quality of Service (QoS) in terms 
of bit-rate and Bit Error Rate (BER), Recently, several ap- 
proaches have been proposed to exploit the MIMO-OFDM 
configuration from a single-user point of view [2] [3]. How- 
ever, one of the main capacities of this structure is that it 
is possible to implement the Space Division Multiple Access 
(SDMA), which means that several users can access the ra- 
dio channel at the same time and at the same frequency. In 
this paper we analyze the design problem corresponding to 
the joint optimization of the receivers and transmitters in 
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a scenario with several parallel communications, assuming 
that Channel State Information (CSI) is available at all the 
terminals. Under this consideration, it is proposed a simpli- 
fied single-user detector for the receivers in Section 2, and a 
joint design for the transmit beamvectors in Section 3. The 
goal is tho minimization of the total required transmit power 
subject to different constraints, such as the maximum BER 
for each user or communication, and/or the maximum power 
for each transmitter. This optimization problem is difficult 
to solve and there are no closed expressions for it. Here we 
propose a heuristic search of the optimal design based on a 
variant of the Simulated Annealing (S A) algorithm and com- 
pare its behavior with a Lagrange-gradient hased technique 
by means of simulation results in Section 4. 

2 SYSTEM AND SIGNAL MODELS 
2.1 System and Signal Models 

We consider the general case in which K users or communi- 
cations are coexisting in the same area, and working in the 
same time-slot and frequency band. The modulation format 
is assumed to be a N-carrier OFDM. In the scenario there 
are several terminals, and each of them Ls able to transmit 
and/or receive, and have more than one Antenna, in this 
paper we address the case corresponding to single-hop point- 
to-point links, in which each communication is assigned to a 
single transmitter and a single receiver, however the exten- 
sion to the case of point-to-multipoint or multipoint- to-point 
(the same communication is carried out by several receivers 
and/or transmitters) links and multi-hop networks is direct. 
In the case we analyze, one terminal can be engaged to more 
than one communication. In the following, the wordy "link", 
"communication" and "user 0 are used as synonyms. 

The snapshot vector signal model for the *th user at the 
nth carrier is [2]: 

y W»»(l) = £ H«*'<*»fcg> 1 jn (t) + oM*»(i) (1) 

jBl 

where r(k) represents the terminal detecting the symbols 
from the kth communication, and *(/) is the element trans- 
mitting the signal corresponding to the fth link. The num- 
ber of components of y£ r and bjf ; is the number of 
an t e nn as of the r(fc)th and t(/)tb terminals, respectively, 
bj* represents the beam vector applied to aSP(l), which is 
the transmitted data in the nth carrier in the rth OFDM 
symbol for the fth communication, where it is considered 
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Figure 1: Typical configuration in a multiuser MIMO- 
OFDM scenario with 3 users or communications. 




Figure 2: General configuration with several terminals 
and p-to-p links. 

that it is normalized: £{|s£ } (t)|*} = 1. £{•} stands 
for the mathematical expectation operator. The matrix 
U<'<i>.>-(*)) represent ^ e MiMO channel for the nth carrier 
between the t(f)th and the r(*)th terminals. Finally, the 
vector ni r(fc)) (l) models the contribution of noise plus inter- 
ferences from outside the system at the r(Jfc)th receiver aL 
the nth frequency, where its associated co-variance matrix is 
*S"*» = E {n™\t)«!,'<'»"(t)}. (■)" stands for complex 
conjugate transpose. 

This signal model can be easily adapted to well-known 
systems such as cellular communications (both the uplink 
and downlink channels, see Fig. 1), and also to the emerg- 
ing reconfigurablo adhoc networks. Pig. 2 shows a general 
configuration with 4 terminals and 5 communications. 

2.2 Single-User Receiver Optimization 
As far as the receiver is concerned , the optimum design 
should be based on a multiuser detector; however, in most 
practical cases this is unaftbrdable since the associated com- 
putational load is too high; therefore in this work we propose 
single-user receivers. The optimum receive beam vector aj 1 * ) 
for the fcth user at the nth carrier is the one that maximizes 
the Signal-to- Noise and Interference Ratio (SNIR) [2]: 

r(»> = ^)) + ^ H («C0.'f*)) b </) b (0" H <«(i).r(*))* (3) 

1-1 

The estimates of the transmitted symbols are calculated as 
follows: #>(*) - dec {.£ fc) "yi r <*>><«)}, where dec{ } stands 



for decision. The SNIR at the detection stage for the Jfcth user 
at the nth carrier is [2]: 

SNTRW = bf> ff Hj'* 1 ' r ^ ,tf Rl i )" l Hi t<i, ' f( * )) bif> (4) 

Taking into account this result, in OFDM it is denned the 
effective BER as the BER av eraged over all the Bubcarri- 
ers: BER<*> = £ Q (yj^SNIRL**), where we have 

assumed that the interferences are approximately Gaussian 
distributed, and depends on the modulation applied to 
each subcarrier (for BPSK, k m = 2). 

3 SIMULATED ANNEALING BASED 
TRANSMITTER OPTIMIZATION 

In last section we have derived low computational cost single- 
user receivers, and now we provide techniques for the design 
of the transmit beamvectors subject to QoS constraints. 
Obviously, a very important parameter is the required power. 
In case of using several transmit antennas! the power used 
for transmitting the information symbol corresponding to 
the nth carrier of the fcth user is proportional to ||bj?|| 2 . 
As far as the QoS requirements is concerned, here we apply 
constraints directly in terms of the effective BER for each 
user, and so we have K constraints. For the fcth user we 
specify which is the maximum effective BER as follows: 

BER< fc >< 7 <*> fc (5) 

Under these constraints, our goal is the mimrnixation of 
the total transmitted power »: Pr ({b^ji^/./JLi) = 

EiLi I^o 1 ffb^ll 2 - This formulation generalizes the re- 
sults presented in [4], where the design was based on a mul- 
tiuser uplink MC-CDMA scenario with one antenna at the 
transmitter and receiver. The QoS constraints were formu- 
lated in terms of the SNIR instead of the BER. In [5] the 
transmit power was stated to be a prefixed value, and the 
goal was the optimization of the mean quality for alt the 
users, so no QoS could be guaranteed for each communica- 
tion, and the physical resources were not used in an optimum 
way as no power allocation was carried out. In [6] these re- 
sults were extended to the most general case and an iterative 
algorithm was proposed, although it was shown that the de- 
duced technique might converge to a suboptxmal solution. 

Besides the previous stated QoS constraints, it could be 
convenient in some RCPnarios, to add constraints correspond- 
ing to the maximum transmit power in some equipments, 
such as the battery-limited MTs in an uplink transmission 
in a cellular system. In this case we generalize this consid- 
eration and define the transmitted power for the sth equip- 
ment as: J><° = £*^ ( ££r 0 ' lib? 1 !!*. ^ our problem we 

add this kind of constraints (not considered in [4], [5] and 
[6]) to some of the traiismitter equipments. Let us define 
X as the set of transmitter elements to which it is applied 
the maximum transmit power constraint. Hence, these are 
formulated as P^ £ P^ax 

There are no closed expressions for the stated constrained 
optimization problem. Some previous works propose gradi- 
ent search algorithms to find a solution [4] or other kind of 
iterative methods [5] [6], but in general, this kind of strate- 
gies may find a local suboptimum minimum due to the non- 
convexity behavior of the optimization problem, as it was 
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shown in [6]. Besides, the gradient techniques may have 
some convergence problems related with the speed and the 
convergence itself [4] [7], and require that the constraints are 
different! able, limiting their application. In this paper we 
propos e the use of the heuristic search algorithm Simulated 
yvnnealiD§ (S A^to find the global opTirnunTpoinT t&Idng'mto 
account- all the constraints previously specified^ 'We suppose 
ffi&TCn^tte solution exists, that is, acofTeclian of transmit 
beamvectors that satisfies ail the constraints simultaneously. 
In case this is not possible, the algorithm will not converge 
to any acceptable design. 

SA is an iterative process able to find the optimum solu- 
tion, even when the problem' tf rior^dfiveiT This"" algoritlim 
~Eas analogies with the annealing of solids in physics as it is 
explained in [8]. In our problem, in each step there is a col- 
lection of beamvectors {b&* , }*fjy;;/ t j_ 1 which is called the 
current solution. Given the current solution (which is equiva- 
lent to a concrete particles arrangement ox state in physics) a 
new solution or collection of beamvectors is proposed. If it is 
"better" than the original one, then it is retained as the cur- 
rent one. On the contrary, if it is "worse", then the proposed 
solution is accepted with a certain probability. This means 
that "worse" solutions may be accepted. This mechanism, 
called "hill-climbing", is important so as to avoid finding a 
suboptimal solution or local mmimiim The parameter that 
controls this acceptance probability is the temperature T t as 
in the case of annealing in physics. The higher the tempera- 
ture, the higher the acceptance probability. The temperature 
is lowered step by step, so that asymptotically only "better*' 
solutions are accepted and a minimum is approached. The 
meaning of "better" and "worse" is related with the defi- 
nition of a cost function /(•) that depends on the transmit 
beamvectors, and which corresponds to the energy of a state 
in the annealing in physics. 

Here we summarize the bask ideas of the SA algorithm 
that we propose to solve the stated optimization problem. 
In this case, we apply the algorithm to a continuous solu- 
tion space. Besides, the cost function /(-) depends on the 
temperature T (6): 

• Cost function definition: 

,(0*.,) = ^'t^ngLf* 

i.r V 'max/ 
where (x) + = x if x > 0 and 0 if x < 0. 

• Proposed solution generation: 

= b?> + w?\ n = 0,... t JV-l 1 Jb = l l ... > A' 
w< fc) ~ Gauss(0,cx£l) (7) 

• Probability of acceptance of the proposed sohdion: 

Prob = exp {-J (/ ({Ef>}) - / ({b<">)) + } (8) 

• System cooling: The temperature is lowered with an 
exponential profile: 

T*-j3T 0 ^ 0.99 (9) 



">**< 1 . 1 . . — 




Figure 3: Performance of the Simulated Annealing in 
scenario 1. 




Figure 4: Performance of the Simulated Annealing in 
scenario 1. Power constraint in MT 1. 



Initially the temperature 7 must be high enough so that 
most of the proposed solutions are accepted. In this paper 
wc run 100 iterations per each value of T. The cunt function 
is equal to the total transmit power plus a quadratic penalty 
term that takes into account if the BERs are greater than 
the required ones, and if the individual transmit powers are 
greater than those specified. As T is lowered, the penalty 
term is increased, and so, we asymptotically avoid solutions 
that do not fulfill the constraints. We make relative com- 
parisons of the BER and transmit powers with the required 
values by means of the log(-) function, as experimentally we 
have observed that it behaves better than absolute compar- 
isons; however, other kinds of penalty functions could have 
been used. The proposed solutions are generated by applying 
independent complex circularly symmetric Gaussian noise to 
the components of the beamvectors. The acceptance ratio is 
monrtorized per each value of T. In case it is lower than 0.1 
for 5 times, then the variance of the Gaussian noise is lowered 
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Figure 5: Performance of the Simulated Annealing in 
scenario 2. 



Figure 6: Performance of the gradient based algorithm 
in scenario 2. 



by means of an exponential rule (of +- 0.95<x?). This helps 
the algorithm iu find with more precision the minimum as 
the temperature is lowered, and so, it increases significantly 
the convergence rate of the technique. 

In the next section we make comparisons between the pro- 
posed SA technique and a gradient solution based on the 
Lagrange multiplier method and the quadratic penalty func- 
tion [4] [7J. In this gradient solution, we only take into ac- 
count the QoS constraints. Under this assumption, it can 
be shown that the minimum transmit power can be achieved 
when the constraints are fulfilled with equality. The tech- 
nique is based on the definition of the Lagrangiau expression 
L Pr + A£f =1 (log 2Mfi£) 3 th e restllts dednced 
in [4], One important problem of this gradient technique is 
that a step-size parameter n must be adjusted, and that, as 
it will be shown In the next section, the speed of conver- 
gence of this technique decreases importantly as the design 
approaches the constraints. Besides, a suboptimal minimum 
may be achieved instead of the global optimum design. 

4 SIMULATION RESULTS AND CONCLU- 
SIONS 

In this section we simulate an uplink channel with 3 MTs and 
1 BS. The OFDM modulation consists of 16 carriers and both 
the MTs and BSs have 5 antennas. The QoS constraints 
in terms of the mean BER are 10" 3 , 10 * 3 and 10" 4 and 
a = 100. In the first scenario, we assume that the path-loss 
is very similar for all the users. In Fig. 3 we show the power 
for the 3 users and the BER, concluding that the algorithm is 
able to find a design fulfilling the constraints when no power 
constraints are applied. The power corresponding to the first 
user is 8.45 W, and the total power is 20.1 W. If we apply 
a power constraint to the first user equal to 8 W, then the 
results are those shown in Fig. 4. We conclude that in this 
case the algorithm allocates 7.6 W to the first user, whereas 
the others increase their corresponding consumption. In this 
Case the global transmit power has increased up to 20.8 W. 

In Fig. 5 and 6 we show the results for an scenario in which 
the third user has a path loss with respect to the first two 



users equal to 12 dB. Fig 5 corresponds to the application 
of SA, whereas Fig. 6 corresponds to the gradient-based 
algorithm with a /x parameter equal to 0.001. We conclude 
that with the same computational load, the SA can fulfill the 
constraints, whereas the gradient based technique decreases 
importantly the speed of convergence as the solution is nearer 
from them. 
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